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Abstract 

We propose the use of bright matter- wave solitons formed from Bose-Einstein condensates with attractive interactions 
to probe and study quantum reflection from a solid surface at normal incidence. We demonstrate that the presence of 
attractive interatomic interactions leads to a number of advantages for the study of quantum reflection. The absence 
of dispersion as the soliton propagates allows precise control of the velocity normal to the surface and for much lower 
velocities to be achieved. Numerical modelling shows that the robust, self-trapped nature of bright solitons leads to a 
clean reflection from the surface, limiting the disruption of the density profile and permitting accurate measurements 
of the reflection probability. 
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1. Introduction 

Solitons are localized self-focusing wavepackets 
that can propagate over long distances without 
change in shape, and emerge from collisions unal- 
tered. Their existence is a common feature of non- 
linear wave equations and consequently they are 
observed in many diverse physical systems includ- 
ing water waves, plasmas, nonlinear optics and par- 
ticle physics pp. Bose-Einstein condensates (BEC) 
in dilute atomic gases are well-described by a non- 
linear Schrodinger equation, known as the Gross- 
Pitaevskii equation (GPE), where the nonlinearity 
results from the atomic interactions (collisions) in 
the gas P]. Such systems can support either dark 
solitons |3|4j or bright solitons |5|6|7j . depending 
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on whether the atomic interactions are repulsive 
or attractive, respectively [5]. Bright matter- wave 
solitons manifest themselves as self-trapped con- 
densates where the usual wavepacket dispersion 
is exactly balanced by the presence of attractive 
atomic interactions. Their self-trapped nature of- 
fers many potential advantages for applications in 
atom optics and interferometry [9] . 

The term quantum reflection refers to the pro- 
cess where a particle reflects from a potential with- 
out reaching a classical turning point and is a direct 
consequence of the wave nature of the particle. Sig- 
nificant reflection occurs when the local wave vector 
of the particle k = \/(^oo 2mU(x)/h 2 ) changes 
by more than k over a distance of 1/fc, where fcoo is 
the wave vector of the particle of mass m far from 
the potential U(x). This requires an abrupt varia- 
tion in the potential U(x), such as can be found in 
the vicinity of a solid surface. The demonstration of 
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quantum reflection from solid surfaces is typically 
performed at grazing incidence in order to reduce 
the wave vector normal to the surface [10I1F . The 
advent of ultracold and quantum degenerate atomic 
samples with large deBroglie wavelengths opens up 
new possibilities to study quantum reflection at nor- 
mal incidence with unprecedented control over the 
atomic motion. Reflection probabilities as high as 
20% have been demonstrated for 23 Na condensates 
incident on a solid silicon surface [H] . More recently, 
the use of patterned silicon surfaces has resulted in 
reflection probabilities of 60% [T3] . 

In this paper we propose the use of bright matter- 
wave solitons to probe and study quantum reflection 
from a solid surface at normal incidence. In section[2] 
we present details of the proposed experimental sce- 
nario and describe how the absence of dispersion en- 
ables precise control of the velocity normal to the 
surface, crucial for the study of quantum reflection. 
The remainder of the paper is devoted to a numeri- 
cal study of the quantum reflection of bright matter- 
wave solitons. In section [3] we present the theoret- 
ical framework for this investigation. In section U 
we highlight the signatures of quantum reflection by 
considering the reflection of a soliton from positive 
and negative step potentials. Finally, in section[5]we 
investigate the quantum reflection from a purely at- 
tractive Casimir-Polder atom-surface potential and 
show that reflection probabilities in excess of 50% 
are feasible from bulk surfaces. 

2. Experimental Scenario 

2.1. Overview 

The proposed experimental scenario for the study 
of quantum reflection is depicted in Fig. [TJ The 
bright matter-wave soliton will be formed in an op- 
tical trap [14115) from an 85 Rb condensate [16] with 
attractive interactions, initially located ~ 5 mm 
from the surface of a highly-polished glass prism. 
Full three dimensional confinement of the initial 
condensate is provided by two intersecting laser 
beams [14]. The first is a focussed 1064 nm beam 
that forms a radial waveguide with typical har- 
monic trap frequencies of up to ~ 100 Hz radially 
and below ~ 5 Hz axially. The second beam, from 
a 50 W 1030 nm laser, is only weakly focussed and 
intersects the first at 90°. This beam provides the 
majority of the axial confinement and has no effect 
on the radial confinement. 




Fig. 1. Schematic of the proposed experimental configuration 
for the study of quantum reflection of bright matter-wave 
solitons. The soliton propagates towards the surface in an 
optical waveguide formed by a focussed 1064 nm laser beam. 
A high power 1030 nm beam intersecting the waveguide at 
90° permits the independent, real time control of the ax- 
ial trapping. An optional repulsive evanescent field can be 
added through the total internal reflection of a 532 nm laser 
beam in the prism. The inset shows the total potential (red) 
experienced by the atoms in the vicinity of the surface as 
the sum of the Casimir-Polder potential (purple) and the 
evanescent field (blue). 



This two-beam configuration is chosen so that the 
harmonic confinement along the waveguide can be 
altered. In particular, by exploiting the deflection 
from an acousto-optic device controlling the second 
beam, the point of intersection of the two beams 
can be altered leading to real time control of the 
trap centre along the waveguide. This in combina- 
tion with the control of the optical power in the 
second beam allows all parameters of the axial har- 
monic potential to be manipulated and permits the 
control of the velocity of the soliton along the waveg- 
uide. For example, following the formation of the 
soliton at the intersection of the two beams, the axial 
trapping beam can then be shifted in position along 
the waveguide causing the soliton to be accelerated 
towards the surface. Removal of the axial trapping 
potential when the soliton passes through the trap 
centre leads to a velocity defined by the displace- 
ment and strength of the axial trapping potential. 
The precise control over the parameters of the axial 
trapping beam afforded by standard optical tech- 
niques yields precise control over the velocity of the 
soliton. 
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2.2. Soliton Formation 

To date, the creation of bright matter- wave soli- 
tons [17] has been demonstrated by three groups at 
ENS in Paris [5], Rice University [6] and, more re- 
cently, JILA [7]. All three experiments use a Fes- 
hbach resonance [18] to switch the atomic interac- 
tions from repulsive to attractive, triggering the col- 
lapse instability of the condensate [19 20 12"I|22| . out 
of which the solitons emerge. The Paris and Rice 
experiments both used a resonance in the 7 Li |F = 
1, mp = 1) state at ~ 700 G, and released the con- 
densates into a quasi- ID optical waveguide in order 
to demonstrate the solitonic character of the attrac- 
tive condensates. 

Here we consider the use of solitons formed from 
85 Rb condensates in the vicinity of the ~ 11 G broad 
Feshbach resonance at ~ 155 G used in the JILA ex- 
periment [7] . The experiment at JILA demonstrated 
the controlled creation of both single and multi- 
ple solitons and revealed interesting soliton dynam- 
ics relevant to the present study. In the case where 
two solitons were created in the radially symmetric 
trap, the solitons were observed to oscillate along 
the weaker axial direction repeatedly colliding in the 
trap centre. The amplitude of the these oscillations 
corresponds to a soliton velocity at the trap centre 
of ~ 0.5 mms -1 . This velocity associated with the 
formation of the solitons for the JILA parameters is 
of the correct order of magnitude required for the 
observation of appreciable quantum reflection from 
a surface [12] . 

2.3. Atom-Surface Potential and Quantum 
Reflection 

At distances greater than the atomic scale, the 
interaction potential of a neutral atom and a solid 
surface takes the form of an attractive power-law po- 
tential. For distances, x, shorter than X/2tt, where 
A is the wavelength of the dominant atomic transi- 
tion, the potential U(x) = — C3/2; 3 , has the form of 
a dipole-dipole interaction between the spontaneous 
atomic dipole and its image [23] . At larger distances 
the potential becomes U(x) = —C4/X 4 due to re- 
tardation of the electrostatic interaction; the well- 
known Casimir-Polder result [24J . 

In the low energy limit (k — > 0), the reflection 
probability R for a particle incident normally on a 
surface tends to unity as R ~ 1 — 2 (34k [25126] , where 
/?4 is the length scale associated with the C4 coeffi- 



cient, C4 = /3|fi, 2 /2m. The observation of quantum 
reflection therefore requires a low incident velocity 
and a weak attractive force to the surface. Such con- 
ditions were first realized using helium or hydrogen 
atoms incident on liquid helium surfaces |27|28] or, 
for solid surfaces, by working at grazing incidence 
[10111] . More recently, quantum reflection at nor- 
mal incidence from a solid surface has been demon- 
strated using a 23 Na condensate incident on a sili- 
con surface at very low velocities [12j . However, in 
the low velocity limit (< 2 mms -1 ) the large spatial 
extent of the (repulsive) condensate leads to anoma- 
lous reflection, as the front reflects before the tail 
producing a standing wave that generates dynami- 
cal excitations that disrupt the atom cloud |29|30] , 
Here we demonstrate that the use of well-localized 
matter-wave solitons will eliminate this problem. 
This coupled to the precise control of the velocity 
of the soliton, stands to take the study of quantum 
reflection to a new level. 



2.4. Evanescent Field 

The experimental configuration also allows for the 
addition of a repulsive (or attractive) evanescent 
field in the vicinity of the surface formed by the 
total internal reflection of a blue (or red) detuned 
laser field within the glass prism. This produces a 
potential which decays exponentially with distance 
from the surface; the decay length being determined 
by the laser wavelength, the refractive index of the 
prism and the angle of incidence of the laser beam. 
When combined with the atom-surface potential, 
the repulsive evanescent field leads to a repulsive 
barrier of finite height, in close proximity to the sur- 
face (see inset in Fig. Q]). Studies of classical reflec- 
tion from such a barrier can be used to probe the 
atom-surface potential [31] . Here the low incident 
velocities of the solitons leads to the peak of the re- 
pulsive barrier being located well into the Casimir- 
Polder regime, unlike in the earlier experiments with 
atoms dropped from a magneto-optical trap [31] . For 
example, we estimate that for a soliton incident with 
a velocity of 10 mms" 1 , setting the barrier height to 
be equal to the kinetic energy of the soliton leads 
to the peak being ~ 170 nm from the surface. More- 
over, the addition of both repulsive and attractive 
evanescent fields can be used to engineer a poten- 
tial in the vicinity of the surface that significantly 
enhances the quantum reflection probability [32] . 
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3. Theoretical framework 

In the limit of ultra-cold temperature the mean- 
field 'wavefunction' of the BEC ifj(r,t) is well- 
described by the Gross-Pitaevskii equation [2] , 
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where m is the atomic mass and the the nonlinear 
coefficient is given by g — Airh 2 a s /m where a s is 
the s-wave scattering length. The external poten- 
tial acting on the system is given by V(r). This 
typically includes a trapping component to con- 
fine the atoms Vt{y). For this we will assume a 
cylindrically-symmetric harmonic trap of the form 
Vt{y) — muj 2 (r 2 + X 2 x 2 )/2 where u r is the radial 
trap frequency and A characterises the axial trap 
strength. Note that the atomic density distribu- 
tion is related to the mean-field wavefunction via 
n(r,t) = |V>M)| 2 . 

Under very tight radial confinement (hio r >> /i 
where /i is the chemical potential), the wavefunc- 
tion can be considered "frozen" in the radial di- 
rection to the non-interacting harmonic oscillator 
ground state. By integrating out this dimension, a 
one-dimensional GPE is obtained in terms of the ax- 
ial dimension x with an effective interaction coeffi- 
cient <7id = g/2nl 2 , where l r = \J h/ mui r is the ra- 
dial harmonic oscillator length. For V(x) = the ID 
GPE supports exact bright soliton solutions given 

by 

#*0 = ^sech(x/O, (2) 

where N is the number of atoms in the soliton and 
£ = 2%/ (m\giT>\N) characterises the soliton width. 

When the potential is weakly-varying in space, 
the ID soliton, to first order, behaves like a classical 
particle. Such classical behaviour has recently been 
studied for matter-wave bright solitons by Martin et 
al. [33], following analogous particle approaches in 
nonlinear optics |34j . Quantum reflection, however, 
dominates in the opposite regime of large potential 
gradients where classical particle analogies are no 
longer valid. The quantum reflection of matter- wave 
solitons from a purely attractive potential well has 
been theoretically studied [35] . In that case a sharp 
switching between the transmission and reflection of 
the soliton was observed, with bound states of the 
well playing a key role in the dynamics. 

In three-dimensions, the presence of radial con- 
finement {ui r > 0) supports a 3D soliton which is 



self-trapped in the axial direction [36] . However, in 
3D a collapse instability exists when N exceeds a 
critical value N c [19]. To avoid collapse effects, we 
will restrict our analysis to N < N c . Note that the 
3D soliton may be prone to collapse during its in- 
teractions with the "surface" . For example, the in- 
teraction of a soliton with an infinite positive wall is 
identical to it colliding with another soliton with a 
7r-phase difference. Such collisions have been shown 
to exhibit regimes of collapse instabilities [37] . 

We simulate the interaction of a bright soliton 
with a surface potential in both ID and 3D using 
the respective GPE. Time-propagation of the GPE 
is performed using the Crank-Nicholson technique 
[38j . In ID the initial soliton solution corresponds 
to the analytic soliton solution in Eq. ([2]), while in 
3D it is obtained with a numerical convergence tech- 
nique whereby the GPE is run in imaginary time 
(i — » —it) from a trial wavefunction 38 . Note that 
when performing 3D simulations we assume cylin- 
drical symmetry such that the relevant coordinate 
space is (x,r). The soliton is initially positioned a 
large distance Aa; to the left of the surface where 
the surface potential is negligible. Since we typically 
consider an axially-untrapped geometry (A = 0), a 
velocity kick to the soliton is required to induce the 
dynamics. This is performed by applying the trans- 
formation ip(r) — » tp( r ) exp{imvx/H), where v is the 
required soliton speed. We typically employ veloc- 
ities in the range v = 0.1 — 5.0 mms -1 . Note that 
when harmonic trapping is present in the axial di- 
rection, as discussed later, a velocity kick is not nec- 
essary as the soliton freely accelerates towards the 
interface. 

Following the JILA soliton experiments, we will 
typically consider a 85 Rb soliton with radial trap 
frequency w r — 2n x 17.5 Hz and axial trapping de- 
fined by A = 0.4. Furthermore, the scattering length 
is a s — —0.6 nm. For such a system, the critical num- 
ber is approximately N c = 2700 [39]. We will there- 
fore employ TV « 2000. 

4. Results 1: Tanh potential 

We will first consider a simplified "surface" poten- 
tial to study the quantum reflection of bright soli- 
tons. We employ a potential step of the form, 
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tanh— 
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By tuning the height Vo and width a of this smooth 
potential, we can control the potential gradient ex- 
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perienced by the incident soliton. Note that this 
soliton-step scenario is analogous to the propagation 
of optical bright solitons at the interface between 
two nonlinear dielectric media of different refractive 
indices [34] . When the soliton is small compared to 
the lengthscale of the step (£ << a), the soliton will 
behave similarly to a classical particle. Since we are 
interested in studying quantum reflection, we will 
consider the opposing regime of sharp tanh poten- 
tials where a < £. 

It is also interesting to mention analogous studies 
of dark matter-wave solitons interacting with step 
potentials. These are localized density dips, sup- 
ported by repulsive atomic interactions, and which 
propagate through an ambient condensate. This va- 
riety of soliton shows no evidence of quantum re- 
flection when incident on a negative step potential 
[4Ti] . Here it is likely that the healing (smoothing) 
of the ambient condensate across the step ensures 
that the dark soliton experiences a smooth effective 
potential, thereby preventing quantum reflection. 



4.1. Negative potential Vq < 

We will first consider the regime Vq < 0, which is 
illustrated schematically in Fig. E^b) (inset) . Clas- 
sically, an incident particle would always transmit 
across such a potential, gaining speed in the process. 

We simulate the impact of the soliton upon the 
tanh potential using the ID GPE. Assuming the 
JILA parameters (outlined above) and N — 1750, 
the soliton has a healing length of £ = 6.4 fim. The 
radial direction (which is assumed but not explic- 
itly calculated) is the harmonic oscillator ground 
state with characteristic width l r « 2.6 /im. We fix 
the height of the potential to be Vb = — 1CP 31 J w 
— 10 ftw r and vary the width a. 

Fig. Ufa) illustrates the dynamics of a soliton in- 
cident on a hard negative step potential [a = 0). 
At low speed [case(i)], the soliton reflects elasti- 
cally from the potential. Since a classical particle 
would transmit across the boundary this is a clear 
sign of quantum reflection. In the opposite regime of 
high speed [case(iv)], the soliton transmits over the 
boundary with negligible perturbation. In between 
these limits [cases (ii) and (iii)] , we observe a regime 
where the soliton splits into a reflected soliton and a 
transmitted soliton at the potential. As the speed is 
increased, we see a gradual reduction of the reflected 
soliton and increase of the transmitted soliton. The 
reflected soliton typically has the same speed as the 
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Fig. 2. Dynamics of a bright soliton incident on a negative 
potential step (Vb < 0). (a) Density plots showing a soli- 
ton incident on a hard (a = 0) negative step for various 
incident speeds (i) v = 0.05, (ii) v = 0.2, (iii) v = 0.6 and 
(iv) v = lmms -1 . (b) Soliton reflection probabilities R as 
a function of incident speed v according to the ID GPE 
(points) and the corresponding non-interacting plane wave 
result (lines). We consider the negative tanh potential of 
Eq. (0 with widths of (from top to bottom) = 0, 0.1, 
0.25 and 1. The inset presents a schematic of the system. 
The step amplitude is fixed throughout to be Vq = — 10 31 J. 



incident soliton, while the transmitted soliton gains 
kinetic energy as it passes down the potential and 
so typically has larger speed. We also note the ap- 
pearance of density fringes as the soliton interacts 
with the step, with the number of fringes increas- 
ing with the incident speed. Similar fringes appear 
in the collisions of two solitons [37] . 

We introduce the reflection probability R = 
Nl/N, where Nl is the number of atoms to the 
left of the potential (obtained via numerical inte- 
gration) following the interaction. In Fig. [2] (b) we 
present the reflection probabilities R as a function 
of incident speed v for various widths a. The results 
of the ID GPE are shown by points. Note that the 
considerable timescale of simulations at very low 
speeds limited our analysis of this regime. 

We have also calculated the corresponding re- 
flection probabilities for a non-interacting plane 
wave. This is performed numerically by dividing the 
smooth potential into many small intervals, each 
of which has constant potential across the interval. 
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The problem is thereby approximated by a series 
of sharp finite potential steps, where the reflection 
at finite step is given by the standard analytic so- 
lution of the Schrodinger equation and the overall 
reflection probability is simply the product of these 
individual reflection probabilities. These plane wave 
predictions, shown by lines in Fig.[5](b), are almost 
identical to the GPE predictions. This contrasts 
strongly with the case of repulsive condensates, 
where, for v < 2mms~ 1 , interatomic interactions 
cause R to fall significantly below the values pre- 
dicted for a single atom [12129] . Over the range of 
speeds presented, the kinetic energy of the soliton 
is sufficiently large that the interactions have neg- 
ligible effect. We do see small deviations appearing 
at low speeds, where the interaction energy of the 
soliton plays a significant role. In this low speed 
limit, the soliton gives slightly enhanced quantum 
reflection over the plane wave solution. 

The quantum reflection curves have identical 
qualitative form, tending towards unity in the limit 
of zero speed, and decaying towards zero as the 
speed is increased. Maximum quantum reflection 
is seen for a sharp step a = 0, as expected since 
the potential gradient experienced by the soliton 
is maximal. Accordingly, for larger widths, the 
potential gradient is reduced and the reflection 
probability decreases. 

An important observation is that the reflection 
probability is very sensitive to the width of the tanh 
potential. For example, at speed v = 0.1 mms -1 , 
R w 0.8 for (7 = 0. For a = 0.25£ this drops to R « 
0.4 and for a = £, the reflection probability is only a 
few percent. Such sensitivity looks promising for the 
use of matter-wave solitons to probe atom-surface 
forces. 

4.2. Positive step V > 

We also consider the soliton dynamics when the 
potential step is positive. Classically, a particle will 
transmit over the barrier when mv 2 /2 > Vq, other- 
wise it reflects. For 85 Rb atoms and Vq — 10~ 31 J, 
the critical speed between reflection/transmission is 
v c w 1.17 nrnis -1 . 

Fig-Elpresents the reflection probability for a soli- 
ton incident on the positive tanh step. The general 
features are that the reflection probability is unity at 
low speed and zero at high speed. In between there 
is a narrow region of transition which is centered on 
the classical critical speed v c . While the GPE results 
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Fig. 3. Reflection probabilities R for a soliton with speed v 
incident on a positive tanh step (Eq. (O) according to the 
ID GPE (points) and non-interacting plane-wave approach 
(lines). We consider step widths of <t/£ = (circles/solid 
line), 0.1 (squares/dashed line) and 1 (stars/dotted line). 
The step amplitude is fixed throughout to be Vb = — 10 31 J. 

vary smoothly from unity to zero, the plane wave 
predictions drop suddenly from unity at v c and de- 
cay smoothly for v > v c . This difference is due to 
the fact that the soliton features a smooth range of 
velocity (momentum) component, unlike the single- 
momentum plane wave. As the step becomes wider, 
the transition region becomes narrower as the soli- 
ton behaves more classically. Recall that for a very 
wide step (a » £) the soliton would indeed behave 
like a classical particle and its reflection curve would 
be a sharp step function centered at v — v c . 

It is important to note that the positive and neg- 
ative steps give very different reflection signatures. 

5. Results 2: Casimir-Polder potential 

We now extend our analysis to a more realistic sce- 
nario in which we simulate the quantum reflection of 
three-dimensional 85 Rb condensates with attractive 
interactions from a plane silicon surface [41j . The 
condensates approach the surface along the common 
(x) axis of circular symmetry. We use two configu- 
rations, a self-trapped condensate with only radial 
confinement which we call "soliton A" and an at- 
tractive condensate propagating in a full harmonic 
trap similar to the bright solitary waves observed in 
the JILA experiment [7] which we call "soliton B" . 

Soliton A has 20 85 Rb atoms and is confined by 
radial harmonic trapping only [u) r = 2n x 17.5 Hz). 
Due to the absence of axial trapping, soliton A is 
only stable due to its attractive inter-atomic interac- 
tions which balance out the dispersion force. Soliton 
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Fig. 4. Schematic diagram of the reflection processes and 
surface potential, (a) and (b) show the initial state for soli- 
tons A and B respectively. Blue ovals represent the solitons 
at their initial position at x = 0. Dotted curves represent 
the longitudinal trapping potentials along the rr-axis and the 
black rectangles mark the surface position at x = Ax. The 
sinusoidal dashed line in (a) represents the initial phase im- 
print. The main features of the Casimir- Polder potential are 
shown in (c) . The thick solid line shows the real part of the 
potential with the dashed line showing the imaginary part. 
The thin vertical solid line shows the position of the surface 
at x = Ax and the vertical dotted line shows the position 
of the potential cutoff. The cutoff is a distance 5 from the 
surface, indicated by arrows. The soliton is represented by 
the hatched area. 

A has longitudinal length ~ 10 /im, radial diameter 
~ 5 /im and an initial peak density uq ~ 10 19 atoms 
m~ 3 [Fig.Sfa)]. Soliton B has 150 85 Rb atoms in a 
three-dimensional harmonic trap (uj x — 2n x 6.8 Hz 
and A = 0.4). Soliton B has longitudinal length ~ 
5 /im, radial diameter ~ 3 /im and an initial peak 
density no ~ 1.4 x 10 19 atoms m~ 3 [Fig. 0(b)]. 

Initially, the condensate is in its ground state, cen- 
tred at x = r = and the silicon surface is posi- 
tioned at a distance Ax away, along the x-axis. At 
time t = the condensates are suddenly propelled 
towards the surface. Due to the different axial ge- 
ometries, we employ different initial conditions de- 
pending on whether we are dealing soliton A or B, 
as shown schematically in Fig. Q] For soliton A, we 
numerically apply a velocity kick to the condensate 
wavefunction at t = [Fig. [4(a)] as outlined earlier. 

For soliton B, we instead employ the harmonic 
trap to set the soliton in motion [Fig. 2(b)], similar 
to the scheme described in section [2] At t = we 
displace the harmonic trap a distance Ax along the 
x-axis, so that the centre of the trap coincides with 
the surface plane. Following the trap displacement 
each atom in soliton B now has potential energy 
V = \m \to 2 (x - Ax) 2 + Lu 2 r 2 ] . The displacement 
causes the soliton to accelerate towards the surface, 
arriving with a velocity v = v x w lu x Ax. The trap 
displacement is typically 5—30 /im, resulting in max- 
imum incident speeds of 0.2 — 1.3 mm s -1 . Since the 
longitudinal radius of soliton B is ~ 5 /im, we cannot 
perform simulations with Ax < 5 /im without in- 
troducing additional complications [42] ■ We do not 
have this problem in the case of soliton A [33] . 
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Fig. 5. Density slices (dark shading indicates high density) 
in x-r plane (axes inset) for 3D 85 Rb condensates with only 
radial confinement reflecting from a planar Si surface at 
incident velocities v = 0.1 mm s — 1 (left, panels (a),(b) and 
(c)) and v = 0.4 mm s -1 (right, panels (d),(e) and (f)). The 
images for the lower velocity case are taken at time t = ms 

(a) , 152 ms (b) and 300 ms (c). The images for the higher 
velocity case are taken at time t = ms (d) , 38 ms (e) and 75 
ms (f). Small arrows indicate the position of density nodes - 
the colour scale in (e) and (f) has been adjusted to highlight 
these. Black rectangles show the position of the surface. 

We now consider the potential due to the Si sur- 
face. For a Rb atom at a distance x' = Ax — x from 
a perfectly conducting Si surface, the interaction 
can be described by the Casimir-Polder potential, 
Vcp(x'), which we model with the single-correction 
function V C p (x') = -C 4 /(x' 3 (:z;' + 3A a /27r 2 )) where 
C 4 = 9.1 x 10~ 56 Jm 4 and A a = 780 nm is the ef- 
fective atomic transition wavelength for rubidium 
[12124129] , This potential has no classical turning 
point and exerts a strong attractive force on the 
85 Rb atoms within a distance of ^2 /im from the sur- 
face. Since Vcp — > — oo as x — > Ax we set V(x) — 
Vcp(S) — i(x — Ax + S)Vi, n for x > (Ax — 5), where 
Vi m — 1-6 x 10~ 26 Jm -1 , and S = 0.15/im is a small 
offset from the surface [29] . The imaginary part of 
the potential models adsorption of Rb atoms by the 
Si surface. A schematic diagram showing the main 
features of the model surface potential is shown in 
Fig. He). 

We now consider the dynamics of the reflection 
process under the combined potential V — Vcp + 
Vt- Fig. [5] shows density slices of soliton A before, 
during and after quantum reflection from the Si sur- 
face at two different velocities. The left panels ((a), 

(b) and (c)) show the reflection process for an inci- 
dent velocity v = 0.1 mm s _1 and the right panels 
((d), (e) and (f)) show the reflection process for an 
incident velocity v — 0.4 mm s _1 . 

First we analyse the lower velocity case. Fig. [5(a) 
shows the equilibrium density profile of soliton A at 
the start of the simulation and its position relative 
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to the surface. We see that the condensate has a 
"cigar" shaped profile. Fig.[5jb) shows soliton A 152 
ms into the simulation. The soliton is just starting 
to move away from the surface having reflected. The 
soliton has been slightly compressed but we see no 
major density modulations because the de Broglie 
wavelength is longer than the soliton's longitudinal 
diameter. Fig. [5](c) shows soliton A after 300 ms. 
The condensate is now considerably more elongated. 
This is because the soliton has lost many atoms (~ 
50%) and hence its axial size, which scales like 1/N, 
grows considerably. The reflection process may also 
excite shape excitations in the outgoing soliton. If 
such excitations are sufficiently large, kinetic effects 
may overcome the attractive interactions and lead 
to the formation of a dispersive wavepacket, as pre- 
dicted by a variational approach [35]. Note that a 
change in the scattering length following the reflec- 
tion could be used to preserve the solitonic nature of 
the condensate. The small modulations observed in 
the density are propagating sound waves. Note that 
the condensate's density profile has no major nodes 
or topological excitations following reflection. This 
is in stark contrast with the low velocity quantum 
reflection of condensates with repulsive interactions 
[29j where lobes and vortices were seen to form. For 
the condensates considered here, the attractive in- 
teractions resist the formation of such features be- 
cause low density regions have higher mean-field en- 
ergies. In this case, the condensate has a low kinetic 
energy (~ 7 x 10 _34 J) compared with the interac- 
tion energy (~ 6 x 10~ 33 J) so interaction effects 
dominate the condensate's behaviour. 

We now analyse the higher velocity case for soli- 
ton A. The soliton starts as before [Fig. Eld)]. At 
38 ms the soliton is just starting to move away from 
the surface [Fig. G2e)] . A density node has formed 
in the cloud (indicated by small arrows) because 
the de Broglie wavelength, ~ 10 /im, is similar to 
the longitudinal diameter of the condensate and 
the condensate's kinetic energy is sufficient to over- 
come the interactions' resistance to sharp features. 
Interestingly, the node is very well defined when 
compared with the non-interacting and repulsive 
interaction cases [29] . It extends all the way to zero 
density and leaves a comparatively large volume 
where we do not expect to find many atoms. This is 
due to the mean field attractive interactions which 
naturally pull atoms away from energetically un- 
favourable lower density areas. Fig. [Sff) shows the 
condensate's density after reflection at t — 75 ms. 
Again we observe an elongation of the cloud because 



the condensate has lost atoms, and hence its ability 
to self-trap. As in the low velocity case, the calcu- 
lations reveal no vortices or topological excitations 
but do show modulations due to propagating sound 
waves. The density node, formed during reflection, 
is still present in the density profile (indicated by 
small arrows). Although the interactions initially 
resist such features, once there is enough kinetic 
energy to form them, this configuration is stable. 

There is a possibility that reflection processes can 
trigger a collapse in attractive condensates if the 
spatial compression at the surface causes the peak 
density to rise above a certain value [36] . Although 
the peak density in these calculations does increase 
by a factor of up to 1.5, we see that this never trig- 
gers collapse for the parameters considered here. By 
analogy to collisions of solitons with ir phase differ- 
ences we would expect collapse instabilities to oc- 
cur for intermediate speeds and for solitons that are 
themselves closer to collapse [37J . 

Fig. [5] shows how the reflection probabilities, R, 
of soliton A (dotted line, square points) and soliton 
B (dashed line, circular points) vary with incident 
velocity v. Both curves show the qualitative trend 
expected for a non-interacting wave packet i.e. R de- 
creases with increasing v. The different propagation 
methods have little influence on the reflection prob- 
ability and the two curves approximately overlap. 
As noted before, the phase imprint method makes 
it straightforward to use lower velocities, hence we 
can reach 0.1 mm s^ 1 without introducing any com- 
plications associated with the initial state touching 
the surface. For this low velocity we predict R = 
0.48. This is comparable with the highest reflection 
probabilities observed in experiments and simula- 
tions with sodium atoms in spite of the fact that the 
rubidium atom is almost four times heavier and no 
nanofabrication technique has been used to enhance 
the reflectivity of the surface [29113] . This is primar- 
ily because lower velocities can be readily attained 
using attractive BECs due to their small spatial ex- 
tent and ability to self-trap. 

We now make a simple comparison between quan- 
tum reflection of condensates with attractive inter- 
actions and condensates with repulsive interactions. 
The inset of Fig. [6] shows R plotted against v for 
simulations of a 23 Na BEC, which is taken to be one- 
dimensional in the rr-direction. The trap frequency is 
lu x = 2tt x 3.3 Hz and the wavefunction is normalised 
such that no = 2.2 x 10 18 atoms m~ 3 as in recent 
experiments [12] . Repulsive condensates do not self- 
trap, hence we use the trap displacement method to 
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Fig. 6. Reflection probabilities, R, plotted against incident 
velocity, v, for 3D 85 Rb condensates reflecting from a planar 
Si surface. The dotted curve with square points is for 85 Rb 
BEC A which only has radial confinement and the dashed 
curve with circular points is for the 85 Rb BEC B with in a 
3D harmonic trap. The inset shows R(v) for a 23 Na BEC 
with repulsive interactions (see text for details). 

propel the BEC along the x-axis towards the sur- 
face. We model the reflection process using the ID 
Gross-Pitaevskii equation with the Casimir-Polder 
potential assuming a s = +2.9 nm and A Q = 590 nm. 

We immediately see that unlike the attractive 
case, reflection probabilities do not continue to rise 
as v — > 0, but rather saturate at a maximum value, 
~ 0.45. This occurs because, at low v, the conden- 
sate spends a long time in contact with the surface 
where the curvature of the wavefunction is high. 
Combined with the repulsive nature of the interac- 
tions, this creates a quantum pressure which drives 
atoms onto the surface. This problem is particularly 
acute for the repulsive condensates we consider here 
because the trap displacements required to attain 
low v are lower than the longitudinal diameter of 
the condensate which therefore starts in contact 
with the surface and remains in contact through- 
out the simulation. An attractive condensate in 
the form of a soliton, on the other hand, prefers to 
avoid contact with the surface. It is, in fact, im- 
possible to create an attractive condensate ground 
state with significant initial surface contact because 
the wavefunction finds sharp features energetically 
unfavourable and subsequently moves away from 
the surface to avoid creating high levels of quantum 
pressure. Attractive condensates have a further ad- 
vantage in that they also move away from absorbing 
potentials - it is possible to reflect an attractive 



condensate from a purely absorbing imaginary po- 
tential which has no variation of its real part. This 
is because the absorbing region of the potential 
draws atoms out of the condensate, thereby raising 
its energy and causing the condensate to move away 
from this region. 

6. Conclusions 

In conclusion we have presented details of a pro- 
posal to use bright matter-wave solitons to probe 
and study quantum reflection from a solid surface at 
normal incidence. We have demonstrated that con- 
densates with attractive interactions have a number 
of advantages in terms of quantum reflection. The 
absence of dispersion as the soliton propagates al- 
lows for precise control of the velocity normal to the 
surface and for much lower velocities to be achieved. 
Our numerical results show that the robust, self- 
trapped nature of bright solitons leads to a clean re- 
flection from the surface, limiting the disruption of 
the density profile. This permits accurate measure- 
ments of the reflection probability which can eas- 
ily exceed 50% and offers a new method to probe 
the atom-surface potential. Such sensitive measure- 
ments in close proximity to a massive object may, in 
the future, present a new method to test the short 
range corrections to gravity due to exotic forces be- 
yond the Standard model [44] . 
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